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Abstract 



The full set of equations governing the evolution of self-gravitating 
spherically symmetric dissipative fluids with anisotropic stresses is de- 
ployed and used to carry out a general study on the behaviour of such 
systems, in the context of general relativity. Emphasis is given to the 
link between the Weyl tensor, the shear tensor, the anisotropy of the 
pressure and the density inhomogeneity. In particular we provide the 
general, necessary and sufficient, condition for the vanishing of the 
spatial gradients of energy density, which in turn suggests a possible 
definition of a gravitational arrow of time. Some solutions are also 
exhibited to illustrate the discussion. 
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1 Introduction 



This work is devoted to the study of dissipative, locally anisotropic, spheri- 
cally symmetric self-gravitating fluids, with particular emphasis in a set of 
physical and geometrical variables which appear to play a fundamental role in 
the evolution of such systems. These variables are the Weyl tensor, the shear 
tensor, the local anisotropy of the presssure and the density inhomogeneity 
The Weyl tensor [1] or some functions of it [2], have been proposed to 
provide a gravitational arrow of time. The rationale behind this idea being 
that tidal forces tend to make the gravitating fluid more inhomogeneous as 
the evolution proceeds, thereby indicating the sense of time. However, some 
works have thrown doubts on this proposal [3]. Further evidence about the 
relevance of the Weyl tensor in the evolution of self-gravitating systems may 
be found in [4] . 

The role of density inhomogeneities in the collapse of dust [5] and in 
particular in the formation of naked singularities [6], has been extensively 
discussed in the literature. 

On the other hand, the assumption of local anisotropy of pressure, which 
seems to be very sensible to describe matter distribution under a variety of 
circumstances, has been proved to be very useful in the study of relativistic 
compact objects (see [7] and references therein). 

A hint pointing to the relevance of the above mentioned three factors in 
the fate of spherical collapse is also provided by the expression of the active 
gravitational mass in terms of those factors [8]. 

Finally, the relevance of the shear tensor in the evolution of self-gravitaing 
systems has been brought out by many authors (see [9] and references therein). 

Now, in the study of self-gravitating compact objects it is usually as- 
sumed that deviations from spherical symmetry are likely to be incidental 
rather than basic features of the process involved (see however the discus- 
sion in [10]). Thus, since the seminal paper by Oppenheimer and Snyder [11], 
most of the work dedicated to the problem of general relativistic gravitational 
collapse, deal with spherically symmetric fluid distribution. Accordingly we 
shall consider spherically symmetric fluid distributions. 

Also, the fluid distribution under consideration will be assumed to be 
dissipative. Indeed, dissipation due to the emission of massless particles 
(photons and/or neutrinos) is a characteristic process in the evolution of 
massive stars. In fact, it seems that the only plausible mechanism to carry 



3 



away the bulk of the binding energy of the collapsing star, leading to a 
neutron star or black hole is neutrino emission [12]. Consequently, in this 
paper, the matter distribution forming the self-gravitating object will be 
described dissipative fluid. 

In the diffusion approximation, it is assumed that the energy flux of 
radiation (as that of thermal conduction) is proportional to the gradient of 
temperature. This assumption is in general very sensible, since the mean free 
path of particles responsible for the propagation of energy in stellar interiors 
is in general very small as compared with the typical length of the object. 
Thus, for a main sequence star as the sun, the mean free path of photons 
at the centre, is of the order of 2 cm. Also, the mean free path of trapped 
neutrinos in compact cores of densities about 10^^ g-cm.^'^ becomes smaller 
than the size of the stellar core [13, 14]. 

Furthermore, the observational data collected from supernova 1987A in- 
dicates that the regime of radiation transport prevailing during the emission 
process, is closer to the diffusion approximation than to the streaming out 
limit [15]. 

However in many other circumstances, the mean free path of particles 
transporting energy may be large enough as to justify the free streaming 
approximation. Therefore we shall include simultaneously both limiting cases 
of radiative transport (diffusion and streaming out), allowing for describing 
a wide range situations. 

It is also worth mentioning that although the most common method of 
solving Einstein's equations is to use commoving coordinates (e.g. [16]), we 
shall use noncomoving coordinates, which implies that the velocity of any 
fluid element (dcflncd with respect to a conveniently chosen set of observers) 
has to be considered as a relevant physical variable ([17]). 

The paper is organized as follows: In the next section we introduce the 
notation and write all relevant equations. Section 3 is devoted to the analysis 
of different special cases. Finally the results are discussed in the last section. 

2 The basic equations 

In this section we shall deploy the relevant equations for describing a dissi- 
pative self-gravitating locally anisotropic fluid. In spite of the fact that not 
all these equations are independent (for example the field equations and the 
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conservation equations (Bianchi identities)) we shall present them all, since 
depending on the problem under consideration, it may be more advantageous 
using one set instead of the other. 

2.1 Einstein equations 

We consider spherically symmetric distributions of collapsing fluid, which 
for sake of completeness we assume to be locally anisotropic, undergoing 
dissipation in the form of heat flow and/or free streaming radiation, bounded 

by a spherical surface S. 

The line element is given in Schwarzschild-like coordinates by 

ds"^ = e'^dt^ - e^dr'^ - (dO^ + sin^edcf)^) , (1) 

where i^{t,r) and A(t,r) are functions of their arguments. We number the 

coordinates: x° — t]x^ — r; — 9; — cj). 

The metric (1) has to satisfy Einstein field equations 

g; = SttT;, (2) 

which in our case read [18]: 




-SttTI = -SttTI = -^(2A + A(A - z>)) 

+^{2u" + ly" - AV + 2^^^) (5) 

-SttTio = --, (6) 

r 

where dots and primes stand for partial differentiation with respect to t and 
r, respectively. 

In order to give physical significance to the components we apply the 
Bondi approach [18]. 
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Thus, following Bondi, let us introduce purely locally Minkowski coordinates 



dr = e'^^'^dt ; dx = e^^'^dr ; dy = rdO ; dz = rsinOdcp. 

Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



7^0 



rp2 
^2 



rp2 . 
^2 ) 



rp3 

^3 



rp3 . 



01- 



Next, we suppose that when viewed by an observer moving relative to these 
coordinates with proper velocity cu in the radial direction, the physical con- 
tent of space consists of an anisotropic fluid of energy density p, radial pres- 
sure Pr, tangential pressure Pi, radial heat flux q and unpolarized radiation 
of energy density e traveling in the radial direction. Thus, when viewed by 
this moving observer the covariant tensor in Minkowski coordinates is 



-q-e 


V 



-q - t 
Pr + e 










p± 





\ 




p± I 



Then a Lorentz transformation readily shows that 

PrU? 2ujq e(l 



7^U 7^U 



Tl = Tl = 



P 



00 



00^ 



rp2 

^2 



rp3 

-^3 



rp2 
^2 



2ujq 

rp3 

^3 ~ 



1 - LO 



UJ] 



(7) 

(8) 
(9) 

Note that the coordinate velocity in the {t,r,9,(f)) system, dr/dt, is related 
to UJ by 



± 



01 



e—T( 



Pr]uje 2 



\ + v 

qe 2 



01 



e^2''e(l -|- uj) 



dr 
UJ — — e 
dt 



{\-v)/2 



(11) 
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Feeding back (7-10) into (3-6), we get the field equations in the form 



1 — uj^ 1 — uj^ 1 — uj Stt \ r I \ 



+ ^ + + (13) 

1 — 1 — 1 — u! Stt \ H r \ 



P, = -^{^(2A + A(A-.)) 
^'^ ^2u" + v'^ - AV + 2^^^) |, (14) 

(15) 
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(p + Pr)u;e 2 2 2\ ^ 2 e(l + u;) A 

1 — oj^ 1 — 1 — UJ QTrr 



The four-velocity vector is defined as 

e 2 uje 2 
«" = (- ^,0,0). (16) 

(1 — U-^ji (1 — 00-^)2 

From which we can calculate the four acceleration a" = u'^pU^ to obtain 

.a, = -«..^ = + 2 ) + + (17) 

For the exterior of the fluid distribution, the spacetime is that of Vaidya, 
given by 

ds^ ^(1- j de + 2dudR - (de'' + sin'9d4>'') , (18) 

where w is a coordinate related to the retarded time, such that u — constant 
is (asymptotically) a null cone open to the future and i? is a null coordinate 

The two coordinate systems (t, r, 6, (f)) and (m, i?, 0, (f)) are related at the 
boundary surface and outside it by 
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(19) 



R = r. 



(20) 



In order to match smoothly the two metrics above on the boundary surface 
= ''"s(^), we require the continuity of the first and the second fundamental 
forms across that surface, yielding (see [19] for details) 



Where, from now on, subscript S indicates that the quantity is evaluated on 
the boundary surface E, and (23) expresses the discontinuity of the radial 
pressure in the presence of heat flow, which is a well known result [20]. 

Eqs. (21), (22), and (23) are the necessary and sufficient conditions for a 
smooth matching of the two metrics (1) and (18) on E. 




(21) 




(22) 



[n]^ = ME- 



(23) 



2.2 Conservation laws {Tj^.^ = 0) 



The energy-momentum tensor (7)- (10) may be written as : 



(24) 



with 



P^ Pr + 2P± 
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p = p + e. 
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9 = 9 + e, 

U = Pr-P± 

and is defined as 



_ A 



= ( r, r, 0, 0) 

with the properties s^m^ = 0, s^s^ — —1, and — qs^. 
We may write for the shear tensor 

with 

1 , V ■ 2ujuj , A , , 2uj' 2uj , 



a — 



ib^nA + z ^) + e-2{uu' + - —)]. 



r 



Then from Tif.^ = 0, using (24), we find: 

p.^u^ + {p + P)e + = Ho.p'J''^ + qal's, 

and 



(p + P)a« + hliq-^u^Sfs + qsp.^u" - P^p + Wf^.^) + (Tapqs'^ + '^Oqs^ = 
Or, contracting (29) with s": 

Pr.^s'' + (P, - P^)sf^ - (p + Px)ss'^ + ^Oq + q-^u" - qs'^s^a^, 

2.3 Ricci identities 

Ricci identitites for the vector Ua read: 

or using 



we have 
1 
2 



1 

■Wa;/? = ttaUf) + aa/3 + -Ohap, 
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2.3.1 Raychaudhuri equation 

Contracting (32) witli -u^ and tlien tlie indices a and /x, we find the Ray- 
chaudhuri equation for the evolution of the expansion: 



q2 

vM" + — + cr«^(j"^ - af^ = -uyRPp = -47r(p + 3P) , (33) 



where 



2.3.2 Constraint equation 

If in (32) we contract first a and /x and then contract with h"'^, we obtain the 
constraint equation expressing a direct relation between expansion 9, shear 
cr"'^ and the heat flux q: 

RppUpN"^ = /i^((7f/ - -ff^^) + (7"^a^ = Stt^s". (34) 

' 3 

2.3.3 Propagation equation of the sheEir 

Contracting (32) with u^h^h'^ we have 

9 9 li" 

(35) 

On the other hand we know that the Riemann tensor may be expressed 
through the Weyl tensor C^^^, the Ricci tensor Rap and the scalar curvature 
R, as: 

p p 1 1 1 

-\%9ocp-^R{5$gan- QaH^I,)- (36) 
Contracting (36) wiih.UpU^h'^h'^ and using Einstein equation (2), we find: 

47r 

Kp^.^P^'^KK = + 47rn^, + y V(P + 3^). (37) 

where E^i, denotes the "electric" part of the Weyl tensor defined by the 
equation (42) below. 
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From (35) and (37) taking into account (33) it follows: 

2 11 

(38) 

2.4 Evolution equations for the Weyl tensor 

According to Kundt and Triimper [21], Bianchi identitites 

RhukS;\ + Rij.uXk;S + RhuSX;k — 0, (39) 

may be written as: 

Cf,uK^;X = RK[fi;u] - -9K[^iR,u]■ (40) 
Then taking into account Einstein equations (2), (40) reads: 



Stt 



Cm^k ;A = 87rT«[^;^] - Y^Kit^T^]- (41) 

In the spherically symmetric case the "magnetic" part of the Weyl tensor 
vanishes (i^a/j = 0), then we have: 

with g^^^ap = gixaQup - gufsgua, ^ixvafi IS the Levi-Civita symbol multiphed by 
and E^'^, the "electric" part of Weyl tensor, may be written as: 

EaP = E{SaSi3 + ^Kji) (43) 



with 



^=£r(U^) (44) 



+ 7. Z— + 



4 2 r 

Contracting (42) with we obtain: 

u'^Cf^.kx = E^kUx - E^xUk-, (45) 
11 



from where it follows that: 

u^'C^J.^^ + <,C^,,^ = + u'^E,,.,o. - u,.,,Ei - u^E^.^,. (46) 
Also, from (31) and (42), we obtain: 

n-xC^J = u^u,aspE'^-afsu^E^,-h^,a''^E^fS+(T.aE;+(T^E:~E^.. (47) 

Replacing (47) into (46), it results: 
A9 



+ apu^E^ + h^^a"^ Eap - a^aE^ - a^aE^ . 
Contracting (48) with h'^ we have: 



(48) 



49 

+ hapcr'^" Ei^„ — (T^aEp — CTk/sE^. 



Kh'pW'C^J.^, = ^E^p-Ufs,,E', + u'^E,,,^h^^h; + 



(49) 



On the other hand 

/i(i/i^«^r,,;^ = (p + p)(a«^ + fM-n^,(< + f/i^) + g,;^/i(i^^ \ (50) 

Feeding back (49) and (50) into (41) we find: 

An 

9Eap + {u'E^^,, - Anu^U^^,, + 47rg«;^)/i^/i^ + YM"p;./ia/3 + Ea^/s = (51) 

■^9 9 
-47r(p + P){aap + -Kp) + 47r(^„a^ + g^Oa) + 47rn^;3« + -/i^) 



Next, contracting (48) with we have: 
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(52) 



The following expressions can also be easily calculated: 



(53) 



Finally, feeding back (52) and (53) into (41) and contracting with /i^ we 
have: 



4:71 



2.5 Weyl tensor, mass function and anisotropy 

For the line element (1) we have: 



^232 ~ 1 ^ 



-A 



2m 



Where the mass function m(r, t) is defined as 



Then from (36), (43) and Einstein equations (2) it follows: 

3777- 

— = 47rp + 47r(Px - Pr) + E 



which in tensorial form reads 

Eap - 4nUaf3 = (— - 4'Kp){SaSf3 + -ha^) 
V o 



(54) 



(55) 



(56) 



(57) 



(58) 
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2.6 Summary 

Equations (28), (29), (33), (34), (38), (52), (55) and (58) are 

+ (p + P)9 + g", = n,^(7°^ + qd's, (59) 



(p + P)a„ + hiiq.^sp + - + H^;^) + a^^qs^ + ^^gs^, = (60) 



e.au"" + ^ + a«;3(7"^ - af, = = -47r(p + 3P) (61) 



= /i;^(af/ - -^;^) + (7°^a^ = ^-Kqs"^ (62) 



2 11 

(63) 



47r 

+ (m'^^^^;^ - ^TTu'^U^^,^ + ATiq^.^)W^hl + yM''p;^/i„^ + ^a«;3 = (64) 

^ Q 

-47r(p + P){(Tal3 + -/ia/?) + ^T^iQadH + ^^Oa) + A'KU^p{a'^ + -h^J 



Ei,xK + «'^aA = -87rg«(<7S + ^/i^) + y(2p + ^P);.K (65) 
-47r(p + P)aa + 47ra''na« - A-ku"" q^.-^hi'^ 
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Eaf3 - 47rn„/3 = (— - 47rp) {SaSp + -hafi) (66) 

In each of equations (59)-(66) there is only one scalar independent com- 
ponent, thus contracting with we may write the equivalent set: 

p* + {p + Pr)e = ^(^ + |)n - - 2qa - (67) 

+ {p + Pr)a + = |g - r - y g (68) 

^* + ?r + T^-a^-«'-^ = -47r(p + 3^) + SttH (69) 
3 6 r 

(f + e)^ = + 12vrg (70) 

£; + 4vrn = -at-a^-^-^ + ^ + ^ (71) 

2 3 r 12 ^ ' 



(47rP, + —){e + -) + {E- AttU + Airp)* = —q (72) 



(£; + Anp - AnUy = — (47rn - E) + Anqi^ + 6) (73) 

r 2 

^=47rp + 47r(Px-Pr)+^ (74) 
with = f aS"", f* = f aU'^ and a" = as", and where the expansion is given 

by 

Then from (27) and (75) it follows at once that: 

(7 ^ 3a;s^ 

- + »=— (76) 
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3 Special cases 



We shall now apply equations (67)-(74) to analyze different particular cases. 

3.1 Geodesic fluids 

If the fluid is geodesic, non-dissipative and locally isotropic, then for bounded 
configurations, it follows at once from (68) and the vanishing of the pressure 
at the boundary, that it should be dust. In this case the vanishing of the Weyl 
tensor implies the shear free condition as it follows from (72). On the other 
hand the shear free condition implies conformally flat as it follows from (71). 
Thus in this special case both conditions are equivalent. For non-geodesic 
fluids this equivalence is not generally true (see below). 

3.2 Locally isotropic perfect fluids 

Let us now consider locally isotropic and non-dissipative fluids {U — q — 
e = 0) and flnd the relations linking the Weyl tensor, the shear and the local 
density inhomogeneity. Although almost all results in this case arc known, 
we think that it is worth while to present them, in order to illustrate the 
general method that will be used later to study more complicated situations. 
Prom (73), we obtain after some rearrengements (with 11 = g' = 0) 

[r^E]^ + r^Anp^ = 0. (77) 

Next, it is convenient to write (72), with the help of (11), (67), (74), and 
(76), as 

[r^E]- + ^[r^E]' = -2Trar\p + P,)^/^^^e'^/^ (yg) 

Implying that the vanishing of the Weyl tensor results in the vanishing 
of spatial gradients of energy density and the shear tensor. 
Let us now assume = 0,then we obtain from (77) 

[r^E]^ = 0. (79) 

Implying, since the Weyl tensor should be regular inside the fluid distribution, 
E = 0. Thus E = and = are equivalent, and either one of them implies 
cr = 0. These results were already known (see [22] and references therein) 
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Next, if cr = it follows from (70) that 

e'' = 0. (80) 

Observe that from the above it follows, using equation (29), that if the 
fluid is conformally flat and satisfies a barotropic equation of state of the 
form Pr — Pr{p), then the fiuid is geodesic {a"' — 0). 

Also, assuming the shear-free condition alone (cr = 0) it follows from (78) 
that the convcctivc derivative of Er'^ vanishes, which in turn means that such 
quantity remains constant for any fluid element along the fluid lines. 



3.3 Locally anisotropic non— dissipative fluids 

We shall now relax the condition of local isotropy of the pressure, and shall 
assume 11 7^ 0. Then from (76) and (73), it follows that 

[r^{E - 47rn)]t + r^4V = 0. (81) 

Implying that the vanishing oi E — 47rn results in the vanishing of spatial 
gradients of energy density. 

On the other hand if we assume the vanishing of p\ then assuming that 
all physical variables are regular within the fluid distribution it follows at 
once that 

E - 47rn = 0. (82) 
Thus E — AttH = and = are equivalent, but neither one of them implies 

(7-0. 

Therefore if we assume the spacetime to be conformally flat {E — 0), then 
the local anisotropy produces inhomogeneity in the energy density according 
to the equation 

(r^n)^ = r^pl (83) 
Next, it follows from (72), with the help of (11), (67), (74) and (76) that 

P(E-47rn)]- +'^[r%E-47rU)]'+87rUuje^''-^^/\^ = -27rar\p+Pr)Vl^e''/^ 

(84) 

implying thereby that the convective derivative of r^{E — AirU) is controlled 
not only by a, but also by 11. 
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If £^ = Airll, then the following link between the shear and the anisotropy 
results 

mue-^ = -o-rVl + P^). (85) 

If the fluid is shear-free and E — Anil then it is either static or locally 
isotropic. Of course in this last case the fluid is also conformally flat. 

3.4 Locally isotropic dissipative fluids in the quasi- 
static evolution 

We shall now relax the condition of non-dissipation by allowing q (for 
simplicity we put e = 0), but assuming that the evolution is slow, which 
means that uj'^ = uj = q = X = i' — and q 0{u) (see [19]). Then from 
(76) and (73), we obtain (in the quasi-static approximation) 

q ir'„-A/2 

+ ^ ^ _4^^t^ (86) 

r 

taking into account that in the quasi-static approximation 

pt^p'e-^ (87) 

we have 

E' + — ^ -Anp'. (88) 

r 

Next, it follows from (72), with the help of (67) 

3E A A / 1 \ 

(89) 

or, equivalent ly 



t'EY = + g (z/' - ^) - ^eVV(p + Pr). (90) 



From (88) it follows at once that conformally fiat and p' = are equivalent 
conditions, which is also true in the perfect (non-dissipating) fluid, in the 
quasi-static approximation. 
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Also, from (89) or (90) it follows that E = does not implies shear-free. 
Indeed, assuming ii^ = in (90) wc have 



q + q[iy 



e^/V(p + P,), 



yielding 



q — re~ 



I 



gA/2+i/ 

2r 



a{p + Pr)dr + p{t) 



(91) 



(92) 



If, we further impose the shear-free condition then, one obtains from (92) 

q = r(3{t)e-'' (93) 

leading to a condition on the temperature, which may be obtained using the 
Landau- Eckart equation 



q^^Kh';(T,-Ta,), 



or 



q=-Ke~{T +^)- 



Using (95) in (93) we obtain 



K Jo 



(94) 
(95) 

(96) 



The two functions P{t) and C{t) are simply related to the total luminosity of 
the sphere and the central temperature, through (93) and (96), respectively. 
A simple model satisfying E — a — wil be next presented. 

3.5 A confer mally flat, shear— free sphere, dissipating 
in the quasi— static regime (with e = 0) 

Prom (5) and = we have 

-A . 1. 1 



SttPx = iu'-X' --] + -, 



(97) 



Then substracting (4) from (97), and considering 11 = 0, we obtain 



,-A 



r^ci + 1, 



(98) 
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where ci {t) is an arbitrary function of time. Substituting (98) into (3) yields 

Snp = -3ci. (99) 
Considering E = with (98) we obtain 

v" + — = 0, (100) 

which has the solution 

e'^/^ = (r^ci + If + C3, (101) 

where C2(t) and c^{t) arc arbitrary functions of t. 
Substituting (98) and (101) into (4) we obtain 

SttP^ = ci [l + 2c2e-(^+'')/^] . (102) 

From (27), condition cr = can be rewriten as 

\ r ) ~ 2r ' 
which after integration becomes 



(103) 



a; = (c4 - l^e^^j re'^/^ (104) 

where C4(t) is an arbitrary function of t. 
From (98), (99), (101), (102) we have 

87r(p + Pr)u = - (2ciC4 - cie^/^) rcge"''. (105) 

Now substituting (98), (101), (105) into (6) we obtain 

87rg = (2C1C3C4 + C1C2) re"*^. (106) 

Next, using (21) and (22) in (98) and (101) we obtain 

2M , , 

ci = (107) 
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and 

C3 = Jl-— (I-C2). (108) 
V '^s 

Also, from the junction condition (23) and from (104) evaluated at the bound- 
ary surface, it follows that 



and 



ci(l + 2C2) = (2C1C3C4 + C1C2) _ (109) 



2M / a;s 2M\ , , 

Solving algebraically the system (108)-(110) for C2, C3 and C4, we can express 
these functions in terms of M, r^, <^s and M. 

We shall further specify our model by assuming 

C2=l, (111) 

C3 = 0, (112) 

implying 

e-^ = e'^ = 1 - (113) 
Using (107), (111), (112) and (113) in (99) and (102) we obtain 

Slip = -SuPr = (114) 
r% 

Then using (111) and (112) in (106) we obtain 

87rg = rcie~'^, (US) 

which is our equation (93) with |^ = /3(t). By virtue of (107) and (109), the 
expression for q becomes 

6Mr (1-^) 

Snq^-^- -r^. (116) 

4 (l-2Mr£ V ; 
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Next, using (107), (109), (113) in (96), we obtain 
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1 2Mr-2 
^ 7^ 



IGTTKrs 



, 2Mr2 
:i )log(l- 



with Tc denoting the temperature at r = 0, and from (104) 
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(118) 



Thus our model represents a conformally flat sphere of fluid evolving 
slowly and shear free, with homogeneous energy density and pressure, satis- 
fying the "inflationary" equation of state p + = 0, with an inward {q < 0) 
heat flux. 



3.6 General case 

First of all, it should be noticed that in the presence of, both, anisotropic 
pressure and dissipation, for Bi warped product spacctimes (which include 
the spherically symmetric case) there exist some restrictions on physical vari- 
ables, provided by the syzygy [23]. However in our case, as it follows from 
the definitions in subsection (2.2) we have | q'^q^, |= {si,q^Y (notice that our 
signature is —2), and therefore the syzygy (see equation (37) in [23]) reduces 
to the identity = 0. 

Then, we obtain from (76) and (73), after some rearrengements, 

[r^{E - 47rn - A'Kqu)]^ + r^Ai^p^ = -47^(gu;)V^ (119) 

or, equivalent ly 

[r^{E - 47rn)]t + 47rrV = A-Kr^q{- + 0), (120) 

if = we recover (81). 

Also, from (72), with the help of (67) we get 

[r3(E - 47rn)]- + ^[r\E - At^U)] + Si^Iiuje^^-^^'h^ 
at 

+47rr'v/r^^e'^/2[g(^ - 2a) - q^] = -2T:ar\p + Pr)Vl^e''''' , (121) 
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which yields (84) in the non-dissipative case. 

From (120) it is clear that the appearance of inhomogeneities in the 
energy-density are controlled by the Weyl tensor, the anisotropy of the pres- 
sure and the dissipation. If we express the heat flux through the temperature, 
the relaxation time and the heat conduction coefficient, by means of some 
transport equation (e.g. [24]) then the aforementioned parameters may be 
related with the creation of such density inhomogeneities. 

Indeed, assuming for example the transport equation [25] 

TKiQl* + Q'' = ^h'^'iTu - Ta.) - -kT' [-^) ,f (122) 

where r, k, and T denote the relaxation time, the thermal conductivity and 
the temperature, respectively. We obtain, putting for simplicity e = 0, and 
using (68) 



<1 = 
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p; + (p + p> + ^(f + e) n 


- K (rt + Ta) 


1 + i 
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(123) 



Then replacing (123) into (120) (with e = 0) we obtain an expression which 
brings out the role of thermodynamic variables in the appearance of density 
inhomogeneities, namely 

[r^{E - 47rn)]^ + A-Kr^p^ 
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p; + (p + p.)a + ^(l + e) 
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- K (Pt + Ta) 


1 + i 


i (2a - 59) + 


K* 2T* ' 
K T 





(124) 

A similar conclusion applies to the shear of the fluid, as it follows from 
(121). 

With the sole purpose of bringing out the role of dissipation in the forma- 
tion of density inhomogeneities, let us present a simple and highly idealized 
model. 



3.7 A dissipative model with E — AttU = 

Since we want to exhibit the role of dissipation in the formation of ihomo- 
geneities we shall assume E — Anil = 0, then from (55) and (58) it follows 
that 
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E - A-kH = ^ m = ^pr^ ^ = I - ^pr^ (125) 
Prom (72), (120) and (76) we obtain 

= q , 126 

r 

and 

3sV2(g + Prco) + {fn-y = 0, (127) 

then combining the Einstein equations (12) and (13) with (125), it follows 
that 

87r(p + pV + 3fi)r 
" 3(1 - f pr^) • ^^^^^ 
We shall further assume the equation of state 

Pr = \p (129) 

then, replacing (129) into (128) we obtain 



:i - f PEr |)^ 

1 - f 

Next, using (127) and (126) it follows that: 



= ^ ' , ■ (130) 

1 - ^pr2 ^ ' 



_'iujq / 1 - |^pr2 



:i31) 



and 



_ -pr ± V(pr)2 - rp'{rp' + 4p)(l - f psr|)2 

(rp' + 4p)(l-fp.r|) ^'"'^ 



In the particular case p = p(i) (but p^ 0), these last equations reduce to 

and ^ 

q = -pre^ (134) 
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From junction conditions it is very simple to express p and p, and thereof 
all physical and metric variables in terms of the total mass M, the velocity 
of the surface the radius and the luminosity of the sphere. Equation 
(131) shows how dissipation produces density inhomogeneity. It is worth 
noticing that both kinds of dissipative processes (e and q) may produce such 
inhomogeneity. 

4 Conclusions 

The study of general relativistic gravitational collapse, which has attracted 
the attention of researchers since the seminal paper by Oppenheimer and 
Snyder [11], is mainly motivated by the fact that the gravitational collapse 
of massive stars represents one of the few observable phenomena, where gen- 
eral relativity is expected to play a relevant role. Ever since that work, much 
was written by researchers trying to provide models of evolving gravitating 
spheres. However this endeavour proved to be difficult and uncertain. Dif- 
ferent kinds of advantages and obstacles appear, depending on the approach 
adopted for the modeling. 

Here, we have established the set of equations governing the structure and 
evolution of self-gravitating spherically symmetric dissipative anisotropic flu- 
ids. For reasons explained in the Introduction, emphasis has been put on the 
role played by the Weyl tensor, the anisotropy of the pressure, dissipation, 
density inhomogeneity and the shear tensor. 

The particular simple relation between the Weyl tensor and density inho- 
mogeneity, for perfect fluids (also valid for locally isotropic, dissipative fluids 
in the quasistatic regime), is at the origin of the Penrose's proposal to provide 
a gravitational arrow of time. However the fact that such relationship is no 
longer valid in the presence of local anisotropy of the pressure and/or dissi- 
pative processes, explains its failure in scenarios where the abovementioned 
factors are present. 

From (120) it is apparent that the production of density inhomogeneities 
is related to a quantity involving all those factors ([r^(£' — 47rn)]t — 47rr'^g(| + 
9)). Alternatively, the appearance of such inhomogeneities may be described 
by means of thermodynamical variables (if a transport equation is assumed), 
as indicated in (124). 

This situation is illustrated in the example provided in the last sec- 
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tion, where density inhomogeneity is produced by dissipative processes alone. 
Thus if following Penrose we adopt the point of view that self-gravitating 
systems evolve in the sense of increasing of density inhomogeneity, then the 
absolute value of the quantity above (or some function of it) should increase, 
providing an alternative definition for an arrow of time. 
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